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We present a theoretical construction for closest-to-classical wave pack- 
ets localized in both angular and radial coordinates and moving on a 
keplerian orbit. The method produces a family of elliptical squeezed 
states for the planar Coulomb problem that minimize appropriate un- 
certainty relations in radial and angular coordinates. The time evolu- 
tion of these states is studied for orbits with different semimajor axes 
and eccentricities. The elliptical squeezed states may be useful for a 
description of the motion of Rydberg wave packets excited by short- 
pulsed lasers in the presence of external fields, which experiments are 
attempting to produce. We outline an extension of the method to in- 
clude certain effects of quantum defects appearing in the alkali-metal 
atoms used in experiments. 
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I. INTRODUCTION 



Coherent states |]J are candidate quantum-mechanical states for probing the in- 
terface between classical and quantum mechanics. Obtaining appropriate coherent 
states for a given quantum situation can be a difficult task, even for apparently simple 
systems. For example, although Schrodinger succeeded in finding exact nonspreading 
coherent states for the harmonic oscillator 0, he was unable to find an analogous 
solution for the Coulomb potential pi. This problem has been discussed by many 
authors, and indeed it is now known for the Coulomb case that there are no exact 
solutions representing localized nonspreading packets following the classical motion 
1, |, |, 0, |, 1,0. 

Given the impossibility of solving the original Schrodinger problem, one can seek 
instead closest-to-classical wave packets for the Coulomb potential without imposing 
the condition that the expectation values follow the exact classical motion for all time. 
It is natural to consider first a simplified version of this problem in which attention 
is restricted to the radial part of the motion, with a fixed (say, p-state) angular wave 
function. Physically, a wave function of this type results when a short laser pulse 
incident on a Rydberg atom excites a coherent superposition of energy states with 
identical angular-momentum quantum numbers. 



It turns out that a kind of squeezed state JTT| called a radial squeezed state (RSS) 
provides a good description of a closest-to-classical radial packet of this type |12|, [13]. 
Such states minimize the uncertainty relation in a set of variables for which the 
radial Coulomb problem takes a form similar to that of an oscillator. An oscillator 
squeezed state is like a coherent state in that it follows the classical motion for all time. 
However, it does not maintain its shape, instead having a position-space spread that 
oscillates periodically |14|]- It can be shown that the corresponding Coulomb RSS has 
radial motion following the closest-to-classical motion along a keplerian ellipse. As 
is characteristic of a squeezed state, the RSS exhibits oscillations in the uncertainty 
product and ratio as a function of time. 

Radial wave packets with p-state angular distributions have been experimentally 
produced using short-pulsed lasers to excite a coherent superposition of Rydberg 
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states [|T^, [H], [T7|, |T^, pt^] . The resulting localized Rydberg wave packets have 
features that match some of those of a classical particle moving in a Coulomb poten- 
tial. For example, the initial motion is periodic with the classical keplerian period. 
While radial Rydberg wave packets initially follow the classical motion, they also ex- 
hibit effects due to quantum-mechanical interference. For instance, at times beyond 
a few classical periods they collapse and undergo a cycle of fractional/full revivals 
2TL E3L UM. Indeed, at still larger times, well beyond the appearance of the first 



full revival, the quantum-interference effects cause radial wave packets to undergo a 
new cycle of fractional/full superrevivals |25j, ^6j. Eventually, at times large compared 
with the superrevival time scale, the wave packets spontaneously decay or lose their 
coherence through other decay processes. 

The RSS description reproduces the primary features in the motion of a Rydberg 
wave packet produced by a short-pulsed laser, including the revival and superrevival 
structures. Moreover, the RSS approach has a number of attractive features. A sim- 
ple example is the shape of an RSS wave packet, which, unlike a simple gaussian, is 
asymmetrical along the radial direction. This asymmetry is desirable in a theoreti- 
cal description of a packet produced with a transform-limited laser pulse exciting a 
superposition of unequally spaced Coulomb energy levels. 

It is natural to ask whether this analytical description of the radial Coulomb 
problem can be extended to incorporate the angular dependence characteristic of a 
classical particle moving in a keplerian ellipse. The issue is of particular interest 
at present because experiments are currently attempting to produce Rydberg wave 
packets that move along elliptical trajectories. 

The purpose of the present paper is to address this question. We provide here 
a planar extension of the RSS construction to a description of minimum-uncertainty 
wave packets moving along a keplerian ellipse, and we discuss some of the key features 
of the solution. Like the RSS, the packets derived here are squeezed states. We call 
them elliptical squeezed states (ESS). 

Unlike purely radial packets, which have p-state angular distributions, experimen- 
tal production of wave packets localized in angular coordinates necessarily involves 
generating a superposition of I states. Superpositions of different I states in a man- 
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ifold of fixed n have been produced [27]. Although a packet of this kind is indeed 
localized in angular variables, the degeneracy of the I states in hydrogen means that 
it is stationary and hence fails to follow the classical motion. 

The production of a wave packet localized in both radial and angular variables re- 
quires creating a superposition of both n and I states. To date, the only experiments 
that have detected periodic motion of a Rydberg wave packet with the classical kep- 
lerian period are those performed with radial packets, which have fixed I. However, 
a technique for producing a superposition of both n and / states has been suggested 
f28fl . The proposal involves making use of a short electric field to convert an angular 
state into a localized Rydberg wave packet. Once produced, such a wave packet would 
move on a circular orbit. If an additional weak electric field is present, a wave packet 
moving on an ellipse of arbitrary eccentricity could be produced. Its motion would 
have features similar to those of the ESS we present here. 

We begin our derivation of the ESS with a discussion of the planar RSS appropriate 
for the keplerian motion in Sec. II. The next task, performed in Sec. IIIA, is to 
obtain appropriate quantum-mechanical operators describing the angular position 
and angular momentum of a particle in a central potential. In the remainder of 
Sec. Ill, we proceed to the construction of a wave packet, called a circular squeezed 
state (CSS), that minimizes the uncertainty product for the angular variables, and 
we discuss some of its properties. 

Section IV provides the construction of the ESS by minimizing the appropriate 
uncertainty relations. At the initial time, the ESS may be written as the products 
of RSS with CSS. Since the ESS naturally involve a superposition over both I and n 
quantum numbers, localization in both the angular and radial coordinates is achieved. 
In this section, we also obtain expectation values for relevant physical quantities and 
discuss the time evolution. The ESS move with classical keplerian periodicity along 
an elliptical orbit of fixed mean eccentricity and semimajor axis. 

Since experiments are often performed in alkali-metal atoms rather than with 
hydrogen, it can be important to allow for modifications due to quantum defects. 
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Some of these are treated in Sec. V. Finally, Sec. VI summarizes our results. 



II. PLANAR RADIAL SQUEEZED STATES 

In this section, we outline the construction of the radial squeezed states for motion 
in a plane. Although the results differ in detail, much of the analysis is similar to ref. 
||13|| , to which the reader is referred for a more complete treatment. 

The RSS construction begins with the classical effective radial hamiltonian for 
the Coulomb potential written in terms of the radial variables r and p r and converts 
it to an oscillator description in terms of new variables R and P. The resulting 
classical problem is then quantized, and wave functions are obtained minimizing the 
quantum uncertainty relation obeyed by R and P. These wave functions are the 
RSS. One motivation for this construction is that the minimum uncertainty of the 
RSS with respect to oscillator coordinates incorporates some of the attractive features 
of oscillator squeezed states, while maintaining the RSS time evolution controlled by 
the Schrodinger equation with a Coulomb potential. 

Classically, the effective radial planar hamiltonian for a particle in a Coulomb 
potential is 

where p r is the radial momentum and E is the energy. For E < 0, the radial motion 
is oscillatory between the values r 12 given by r 12 = a(l ±e), where a = 1/2\E\ is the 
semimajor axis of the orbital ellipse and e = (1 — 2l 2 \E\) l l 2 is the eccentricity, with 
I the constant classical angular momentum. The classical orbital period is given by 
T cl = -n/{2\E\fl\ 

The conversion to an oscillator description is via the variable change 

fl=--~, P = Pr ■ (2) 

r r 

In terms of R and P, the equation H = E becomes 

\P 2 + \l 2 R 2 = \e 2 , (3) 
which has the form of an energy equation for an oscillator. 
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Promoting R and P to the status of quantum-mechanical operators R = - — ^ 
and P = p r = —i(d r + l/2r), where [r,p r ] = z, we obtain the commutation relation 
[-R, -P] = —i/r 2 . The RSS are the wave functions ip{r) minimizing the corresponding 
uncertainty relation 

ARAP > (4) 

at fixed time. 

The minimum-uncertainty wave packets are given by 

i/j(r) = N\r a exp[— 7or — ijir] , (5) 

where the normalization constant is Ni = (27 ) a+1 /[r(2a + 2)] 1//2 . The parameters 
a, 70, and 71 are given in terms of the squeezing 

AR 2(AR) 2 
and expectations of the operators 1/r and p r by 



a = \~\ > 7o = \(\) > 7i = -<Pr> ■ (7) 



The wave functions ^(r) form a three-parameter family of planar RSS. The radial 
part of the classical motion is uniquely determined by two parameters specifying 
a point in the r-p r phase space. However, three parameters are needed to fix the 
quantum-mechanical solution ip(r) because, in addition to specifying (r) and (p r ), it 
is necessary to provide information about the spread of the packet. This can be done, 
for instance, by fixing the mean energy (H). 



The RSS discussed in refs. [0, |T3| can be used to describe a radial Rydberg wave 
packet with p-state angular distribution produced by a single short laser pulse. The 
corresponding three RSS parameters for this situation are determined by matching 
the RSS expectation values for (H), (r), and (p r ) to mean values for the Rydberg wave 
packet at the outer apsidal point. Three experimental quantities control this match. 
The frequency of the laser pulse determines the mean energy of the wave packet. 
The time delay after the excitation of the wave packet determines the location of the 
packet. Finally, the duration of the laser pulse fixes the spread of the wave packet. 
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Our purpose in the present work is to obtain wave packets (the ESS) that follow 
the classical keplerian motion along an ellipse. This implies incorporating nontrivial 
angular dependence into the analysis, so quantities such as the mean energy depend 
on both radial and angular coordinates. Again, certain expectation values specify 
the ESS parameters. However, the coordinate dependence is sufficiently intertwined 
that it is impossible to treat independently the determination of the three planar 
RSS parameters a, 70, and 71 of Eq. (|5|). For this reason, we defer until Sec. IV a 
discussion of this topic. 

III. CIRCULAR SQUEEZED STATES 

The goal of this section is to obtain the circular squeezed states, which can be 
viewed as squeezed states for angular coordinates in the plane. Sec. IIIA contains 
some remarks about the choice of angular coordinates. Sec. IIIB derives the CSS, 
while Sec. IIIC examines properties of the CSS. 



A. Quantum Angular Variables 

The effective angular hamiltonian for planar motion is the square L 2 of the angular 
momentum operator L. There are, however, subtleties involved in the choice of the 
corresponding angular coordinates to use in the quantum description. This subsection 
briefly discusses some of the issues directly relevant to the construction of the CSS 
and provides the definitions and conventions we have adopted. For more details and 



a guide to the early literature on this topic, we refer the reader to the review p9 |. 

Among the possible choices for the angular coordinate are the continuous variable 
cf) c with —00 < <fi c < 00 and the periodic variable <p p with —n < (f) p < ir. Both 
these choices present problems at the quantum level. The continuous variable C is 
not periodic, and hence is not an observable in the Hilbert subspace of functions for 
which L = —idcj, is hermitian under a conventional inner product. In contrast, 4> p is 
periodic but is discontinous at (f) p = tt, which causes unusual effects with derivative 
operators such as the angular momentum. For instance, direct substitution into the 
commutation relation [<f> p , L] generates 5-function contributions. 
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These difficulties can be avoided by choosing instead angular coordinates that 
are both periodic and continuous. However, a single such quantity is insufficient to 
specify uniquely a point on a circle, since periodicity implies the existence of extrema 
and hence no single quantity can be one-to-one. One possibility |30]| is to use two 
angular coordinates, cos0 and sin0, where is either C or P . Classically, this is 
a natural choice since it corresponds to the identification (x,y) — > (cos 0, sin 0) on a 
unit circle. 

The corresponding quantum operators can be defined through their matrix ele- 
ments in the Hilbert space. We take as the inner product the definition 

(v>iiv>2>= r d<f>r^2 . (8) 

J —IT 

Matrix elements of the quantum operators cos and sin are then given by insertion 
of the coordinate functions cos0 and sin0 in the inner product (|8|). 

- — ■ 2 - 2 

The angular operators obey the nonlinear relation sin0 + cos0 = J, where / is 
the identity operator. Note that using one of the two angular coordinates along with 
the sign of the other suffices to determine a unique location on a circle. However, this 
choice is not smooth. Taken together, the pair of operators cos 0, sin have all the 
features needed for a quantum- mechanical treatment. 

The quantum angular-coordinate operators cos0 and sin0 and the angular mo- 
mentum L are intertwined by the commutation relations 

[cos0, L] = — i sin0 , [sin 0, L) = i cos , [cos 0, sin 0] = . (9) 

The uncertainty relations following from these equations are 

Acos0 AL > ||(srn0}| , (10) 

Asm0AL > ||(cos0}| , (11) 

and 

Asm0Acos0>O . (12) 

The last of these relations shows that cos and sin can be specified simultaneously 
to arbitrary precision, which is intuitively reasonable since both quantities are needed 
to determine the location of the quantum particle. 
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For simplicity, in the remainder of this paper we adopt the usual convention of 
writing cos as cos and sin as sin 0, leaving the context to determine whether an 
operator or a function is intended. 

B. Construction of CSS 

We begin by seeking a minimum- uncertainty state that is centered about = 0. 
Circular symmetry then implies the conditions 

(sin0)=O , (cos0)>O . (13) 

This means that the uncertainty relation (|ll|) provides the only nontrivial constraint. 
A wave packet minimizing this relation must obey the differential equation 

(sin0) X = -^(L-(L»x , (14) 

where the squeezing 6 in the angular coordinates is a real constant given by 

1 _ Asin0 2(Asin<^ 2 



5 AL | (cos 0)| 

The solution to Eq. fll4|) is 



> . (15) 



X(0) = iV 2 exp(5cos0 + i/30) , (16) 

where we have defined the real parameter (3 = (L). Note that (3 must be integer for 
the x(0) to be single valued. This is a special case of the general result that minimum- 



uncertainty angular packets must have integer angular-momentum expectations [pi 
The normalization constant N 2 is given by 



*isw • (17) 

where I n {z) is a modified Bessel function of the first kind. The two-parameter family 
of states ( [IB] ) are the CSS. They have previously entered the literature in the context 
of uncertainty relations for phase and angle variables |29[]. In addition to (L) = f3 
and (sin#) = 0, already imposed above, we find 

^=W)-° • (18) 
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as expected. 

By construction, the CSS have magnitude peaking at = 0. It is shown in 
Sec. IV that this is appropriate for use in deriving packets moving along an ellipse 
aligned with outer apsidal point at = 0, which is the orientation most likely to be 
relevant in future experiments. However, modified CSS with magnitude peaking at 
some angle = 0o instead would be relevant for packets moving along an ellipse of 
arbitrary orientation. The remainder of this subsection concerns these modified CSS. 

We first remark that imposing (cos0) < instead of Eq. fll3|) results in a sign 
change of the inverse squeezing 8 and hence in a packet of the same shape as Eq. 
flOp but centered about = ir. Another possibility is to seek a configuration with 
(cos0) = instead of fll3|) , whereupon the interesting constraint becomes Eq. ( pTT|) 
instead of fllPl). The solution is again of the form ([ll]) but with sin replaced by cos 
in the exponential, representing a packet centered at = ±7r/2 according to the sign 
of (sin0). 

One approach to the construction of a packet centered at an arbitrary angle 0o is 
to take advantage of rotational covariance. Implicit in the use of cos0 and sin0 as 
angular variables is a choice of origin for 0, or, equivalently, a choice of orientation 
for the coordinate axes in the plane. However, the effective angular hamiltonian L 2 
has U(l)-rotational invariance corresponding to independence of the choice of origin. 
A rotation by an angle 0o may be implemented by a translation in the continuous 
coordinate, C — > C — 0o- This leaves invariant all physical observables since the 
inner product @ is invariant on functions of the periodic coordinates sin# and cos 6. 

The conditions fll3D and the solution (16) break this U(l) symmetry. A U(l) 
rotation of (|T6|) by — > — 0o produces a different packet with maximum centered 
at 0o : 

X(0, 0o) = N 2 exp [5 cos(0 - O ) + i(3{<f) - O )] , (19) 

Instead of Eq. (|T0D, this packet minimizes an uncertainty relation between a new 
angular-coordinate operator cos(0 — 0o) and the angular-momentum operator L. 

An interesting issue is whether there is a meaningful way to minimize both (10) 
and (0). Simultaneous extremization of both relations is impossible. Furthermore, 
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various combinations could be chosen for minimization, producing a variety of inter- 
polations between the solutions found at = 0, ±|vr, tt. However, there is a unique 
prescription for treating both relations together such that the packet produced is the 
rotated version fll9|) of (|T6|). For any given packet, define the quantities 

/i 2 = [Acos(0-0 o )] 2 (AL) 2 

- i|(sin(0-0 o ))| 2 , 
v 2 = [Asin(0-0 o )] 2 (AL) 2 

- I|(cos(0-0 o ))| 2 • (20) 



Requiring that the quadratic combination Q = fi 2 + v 2 be the constant Qo given by 
the case 4>o = yields the solution (|19D . The value of Qo is 



o - P 



7 (25), 



<5 2 \I (25) + h(25) 



(21) 



C. Features of CSS 

The form ( |16|) of the CSS can be used to obtain second-order expectations. We 
find 

<«-'*>= S • (22) 



and 



c^S+^N^ ' (24) 



Note that (sin 2 0) + (cos 2 0) = 1, as required. In the limit 5 — > 0, (sin 2 0) and (cos 2 ■ 
both approach |, while as 5 — > oo, (sin 2 <f>) — > and (cos 2 0) — > 1. 
The uncertainties in the angular coordinates and momentum are 



>2 /„(2^)/;(2i) - [/ 1 (2a>] 2 



w) - W ' (25) 



,2 1 h(2S) 1 
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<-MiH«^> ■ (27) 

These confirm that the uncertainty product ( |TTD is minimized. 

The uncertainty relations reveal that the CSS parameter S represents the angular- 
momentum spread of the solution x{4>)- As 5 decreases so does AL, and hence 
the angular-coordinate width increases. In the limit 5^0, the normalization con- 
stant N2 — > 1/V27T, so tends to an angular-momentum eigenstate Yp((j)) = 
(27r)~ 1//2 exp(i/30). This state has AL — > 0, consistent with the limit 5 — > of Eq. 
( p7|) . Moreover, (sin0) — > also, so both sides of the uncertainty relation (|10|) vanish 
in this limit. This avoids the appearance of infinite uncertainty in Acos0, which is 
impossible for a bounded function. Similarly, as 5 — > both sides of Eq. (|TT|) vanish. 

The interpretation of 5 as the angular-momentum spread is confirmed by an ex- 
pansion of the packet ([16]) in the vicinity of = 0: |x(0)| 2 ^ exp[25(l — + • • •)]. 
This shows that to leading order |x(^)| 2 has a gaussian dependence proportional to 
exp(— 5(f) 2 ), with angular-coordinate standard deviation a = l/y/25. It also implies 
that a CSS of finite angular-coordinate width satisfies | (cos 0) | < 1 along with the 
conditions (|T3|), which differs from the corresponding classical particle located on the 
unit circle at cos</> = 1, sin0 = 0. 

Given the initial angular location, a CSS is specified by two quantities 5 and (3. 
Compared to the corresponding classical problem, an extra parameter is needed to fix 
the quantum solution. Classically, the motion on a circle is determined by the initial 
values of sin0, the sign of cos0, and L. At the quantum level, the initial angular 
position is specified by (sin0) and the sign of (cos0), while (3 = (L) establishes 
the angular momentum. The extra parameter 5 controls the angular spread of the 



packet and can be fixed through Eq. (24) by giving AL. Evidently, in an experiment 
attempting to produce a wave packet localized in the angle 0, it is insufficient to excite 
a state of definite I. A single laser pulse alone therefore cannot excite an angular wave 
packet. A superposition of I states with a spread AL must be produced by turning 
on additional fields that excite more than one value of 

Plots of some normalized CSS as a function of </> are presented in Fig. 1. The 
configurations all have (sin0) = 0, (cos0) > 0, and j3 = (L) = 30. Three different 
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examples are shown, with AL = 0.5, 1.5, and 2.5 corresponding to 5 ~ 0.804, 4.757, 
and 12.753, respectively. As expected, increasing S decreases the angular-coordinate 
spread. 

IV. ELLIPTICAL SQUEEZED STATES 

This section discusses the elliptical squeezed states. Sec. IVA constructs them and 
explains how a given ESS is specified in terms of physical quantities. The issue of 
time evolution is examined in Sec. IVB, where we show that the ESS move along an 
orbit close to a classical ellipse. In Sec. IVC, features of the ESS solutions and their 
relationship to other approaches are considered. 

A. Construction and Specification of ESS 

The CSS solution ( |I"6D is a function only of the angular coordinates and has neither 
time dependence nor dependence on the energy quantum number n. One solution of 
the full planar Coulomb problem is the the product of a CSS with a radial energy 
eigenstate of given n. However, this produces a stationary state. 

The desired closest-to-classical packet moving on a keplerian orbit can be obtained 
by taking advantage of the separability of the full hamiltonian. Note that the un- 
certainty relation (f|) is independent of I because the uncertainty AR involves the 
combination R — (R). This implies it is possible to minimize simultaneously the un- 
certainty relations @ and (|T0|) at a given time. We can therefore take as an initial 
state the product of an RSS 4>{r) and a CSS x(0), giving 

i 

= Uw(2 a + 2) J r"exp[-( To + iTl )r + *W + W ] ,(28) 

where we have substituted from Eqs. @ and flUf). This is a normalized five-parameter 
family of ESS. They represent minimum-uncertainty packets localized in both radial 
and angular coordinates and evolving in time. As discussed in Sec. Ill, the choice 
of the initial angular-coordinate location is implicit in the CSS construction and is 
taken to be (sin0) = and (cos0) > 0. 
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Expectation values of operators in the ESS can be calculated analytically using 
Eq. fl2"8|). Some physically useful expectations are as follows: 

<r> = — , (r , )= (« + l)(2° + 3) (2Q) 

7o ^ 27 



(Pr> = -7i , ^ = ^ + ^ ' ( 3 °) 



(sin0) = O , (cos0) = ^|>O , (31) 

Jo (2d) 

W = /5 , < L2 > = fr^ + /? 2 . ( 32 ) 

{H) ~ 2(2a + l) + T + a(2a + l) fe^j " (33) 
The expectation value for the energy (if) is obtained using the planar Coulomb 
hamiltonian and depends on all five of the parameters associated with the component 
RSS and CSS. 

The associated uncertainty products are given by 



ArAp, = W^t- (34) 



AMAL = i§S (35) 



1 / (2d)/ 1 (2 d)/((2d) -[/ 1 (2d)] 3 
[/o(25)] E 



a cos 9al = \ oy i n :: i ;z\ [ ■ w 



The ESS is not a minimum-uncertainty state in r and p r , which is as expected since 
it is constructed to minimize the uncertainty relation (||) instead. For large values of 
a, however, ArAp r — > |. 

We consider an initial configuration where the wave packet is located at the outer 
apsidal point of an elliptical orbit. This choice is consistent with the experimental 
situation, where the uncertainty product for radial Rydberg packets reaches a mini- 
mum near the outer apsidal point |^TJ . The constraints imposed on (sin 0) and (cos (f>) 
ensure that the ellipse has semimajor axis aligned along the x axis of the coordinate 
system. The parameters a, (3, 70, 71 and 6 can be fixed in terms of the spread AL 
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in the angular momentum and expectations of the radial coordinates (r) and (p r ), of 
the angular momentum (L), and of the energy (H). 

Classically, (r), (p r ), and (L) determine the mean location and initial velocity of 
the packet. At the quantum level, two additional physical conditions must be given 
to determine completely the ESS. These conditions represent the initial width of the 
packet in the radial and angular coordinates. For the first, we set the energy expec- 
tation equal to the mean energy of a Rydberg packet consisting of a superposition of 
n states centered on the value n. A packet of this type can be produced by excitation 
with a short laser pulse tuned to the mean energy E n = —\/2{n — 1/2) 2 . If the 
excitation occurs in the presence of external fields, a superposition of / states can be 
achieved as well. The second extra condition is the spread AL in this superposition. 
The precise field arrangement will determine AL, although in practice it may be 
difficult to determine a priori for a given experiment. 

Denote by / the average value of I for the superposition of I states. Then, we 
choose to set the expectation of the coordinate r equal to a radial distance to the 
outer apsidal point written in terms of these average quantum numbers, 



r ou t = {n-\Y\l + 



• (3-) 



\ (n- 



The expectation of the radial momentum is set to zero. These choices are not uniquely 
enforced by experimental considerations but do ensure that both RSS and CSS are 
readily recovered in suitable limits of the ESS. 

The full set of conditions sufficient to fix the five ESS parameters is therefore as 
follows: 

(r) = r out , (p r ) = , (L) = I , 

(H) = E n , y[{p) - (L>2 = AL . (38) 

These determine a, ft, 70, 71, and 5 in terms of the three quantities n, I, and AL 
whose values depend on a particular excitation scheme using a short-pulsed laser in 
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the presence of external fields. 



B. Evolution of ESS 

With parameters matching a Rydberg wave packet at the outer apsidal point 
on an elliptical orbit, the ESS (pgj) may be taken as a minimum- uncertainty initial 
solution of the time-dependent Schrodinger equation. Since (p r ) = and (L) = [3 
by construction, for /3 > the wave packet will begin to move in the direction of 
increasing <j). The geometry of the ensuing orbit depends on the values of n, I, and 
AL. In the limit I — > n — 1, the orbits for the ESS should become more circular in 
shape, and the wave packet should propagate at a fixed mean radial distance from the 
nucleus. As I — > 1, the orbit should become highly elliptical, and the inner apsidal 
point should approach the origin. In this case, radial oscillations should occur as the 
particle passes close to the nucleus, as is observed for Z = 1 RSS fll2| . 

The time evolution of the ESS may be studied by expanding \l/(r, <f),t) in eigen- 
states R n i(r) of the energy E n = — l/2[n — 1/2) 2 and Yi((p) of the angular momentum, 
given by 



R nl (r) = Nr^e-^Lf\ hl (^-j 



W) = iexpH/0) , (39) 

V Z7T 



where N is a normalization constant, n is a positive integer, and I is a positive or 
negative integer satisfying \l\ < (n — 1). The expansion of the ESS is 

^(r,0,t) =J2cniRm(r)Y l (<j ) )e' lE " t . (40) 

n,l 

Note that the time dependence of the phase means that the ESS is separable only at 
t = 0. 

The expansion coefficients c n \ = (\l/(r, 0, 0)|i? n ;(r)Yi(0)) may be calculated using 

Eq. ( ^8|) as the initial wave function \l/(r, 0, 0). The result for / > is 

1 r i 

2 T , / r, \ 21+2 



Cnl 



(27o 



,2a+2 



I (25)T(2a + 2) 



T(n-l) 
(2n - 1) \n-\) T(n + l) 
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X h U-lJ mp + 2l + inn-l- P ) \n-\ +1 ° + l11 ) 
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As an example, consider a wave packet with n = 45, I = 30, and AL = 2.5. Using 
these values and Eq. (|38|), we obtain the ESS-parameter values a ~ 57.408, (3 = 30, 
7o ~ 0.01697, 71 = 0, and 5 ~ 12.753. The series in (fH]) may be well approximated 
by truncating the sum to a finite number of terms centered on n and I. This permits 
the ESS to be plotted as a function of r and <ft at different times t. 

Figure 2a shows the initial ESS defined at t = 0. It is located on the positive x 
axis at the outer apsidal point and is evidently localized in both radial and angular 
coordinates. The radial distance to the outer apsidal point is r out ~ 3443 in atomic 
units. 

The classical orbital period for motion on a keplerian ellipse is T c \ = 27r(n — 1/2) 3 . 
With n = 45, we obtain T d ~ 13.4 psec. Figure 2b shows the ESS at the time t = |T cl . 
It has moved in the direction of positive <fi and is spreading along the elliptical orbit. 
Since the packet moves slowly near the outer apsidal point but faster near the inner 
apsidal point, in accordance with Kepler's laws, it has traversed less than one third 



of the orbit at t — \T C \. 



Figure 2c shows the ESS at t = \T cl . It is moving rapidly and has substantial 
spread along the elliptical orbit. Since Z = 30 is relatively large for a state with 
n = 45, the ESS remains at some distance from the nucleus and therefore exhibits 



none of the radial oscillations seen for RSS p states |12|]. The radial distance to the 
inner apsidal point is r- m ~ 517 a.u. 

Figures 2d and 2e show the wave packet at the times t = |T c i and t = T c \, 
respectively. The motion is slower again as the outer apsidal point is approached, 
and the packet becomes more localized again. At t = T c \, the wave packet closely 
resembles the initial wave packet. However, the motion is not exactly periodic. For 
times t 3> T c i, the wave packet collapses and a cycle of revivals and superrevivals 
commences. 

The localization in the angle <fi oscillates as the wave packet goes through a classical 
orbital cycle. At t — 0, it is highly localized in angle, while at t = hT c i, the wave 
packet is extended around much of the orbit. By t = T c \, the wave packet has localized 
again. This oscillation is characteristic of a squeezed state. The radial uncertainties 
are oscillating also, although not as noticeably as for a p-state radial wave function. 
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In Fig. 3, two ESS with different values of I and hence different orbital eccentricities 
are compared. Both packets have n = 45 and AL = 2.5. They are viewed at t — ~T c i 
from a point on the positive x axis looking towards the nucleus. Figure 3a, with 
/ = 30, is the same as Fig. 2c but viewed from a different perspective. The ellipticity 
of the orbit is manifest. Figure 3b shows a different ESS with / = 40, for which 
a ~ 20.412, f3 = 40, 7 ~ 0.00752, 7l = 0, and 5 ~ 12.753. Since AL is unchanged, S 
remains the same. As expected, the orbit is closer to circular. For this case, r in ~ 1113 
a.u., which is more than twice the value for the / = 30 wave packet in Fig. 3a. 

C. Features of ESS 

The literature contains examples of quantum wave packets other than the ESS 
that are nonetheless localized on classical elliptical orbits. In this subsection, we 
discuss a few additional features of the ESS and provide a comparison with some of 
these other approaches || || . 

The methods of refs. [§, |9j both involve the Runge-Lenz vector. In the classical 
Coulomb problem, this vector is a conserved quantity with magnitude equal to the 
eccentricity of the orbit and with orientation along the semimajor axis pointing from 
the focus to the inner apsidal point. At the quantum level, the Runge-Lenz operator 
in atomic units takes the form 

A=\{pxL-Lxp)- r - . (42) 

In ref. ||, the Runge-Lenz operator is scaled by a factor involving the energy oper- 
ator at a fixed value of n. Within a fixed n manifold, the Coulomb problem exhibits 
an SO (4) symmetry. Coherent states for the SO (4) symmetry can be constructed 
using the scaled Runge-Lenz operator and the angular-momentum operator. The re- 
sulting wave packets are localized on elliptical orbits. However, they are stationary 
states that do not follow the classical motion. 

A similar approach is taken in ref. ||. Planar motion is considered, and a subset 
of the operators is used to form an o(3) algebra. The operators A x , A y , and L satisfy 

[A x ,A y ] = -2iHL , (43) 
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where H is the hamiltonian. The associated uncertainty relation is 

AA X AA y > | (HL) | . (44) 

Quantum solutions minimizing this relation are found. The resulting states are also 
stationary and are elliptical in shape. By taking a gaussian-weighted superposition 
of these states with different values of n, a localized wave packet is produced that 
moves along a keplerian orbit and follows the classical motion. In the limit of large 
angular quantum numbers, Eq. ( fPj) approximately holds for these superpositions. 

The ESS are different from the solutions of refs. || |9|. They naturally contain 
a superposition of n and I states, and they are minimum- uncertainty solutions for 
both radial and angular operators. The shape of the elliptical orbit is determined by 
the parameters of the ESS, which in turn depend on the experimentally determined 
parameters n, I, and AL. 

A natural question to ask is whether the ESS also provide a minimum-uncertainty 
solution of Eq. ([H]). For the example plotted in Fig. 2, the values of the ESS param- 
eters given in Sec. IV can be used to determine numerically the left and right hand 
sides of (||). We find AA X AA y ~ 0.1214, while \(HL)\ ~ 0.0099 in atomic units. 
Defining the quantity 

z - \m\ ■ ( ] 

we find Z ~ 11.26 for the ESS wave function in Fig. 2a. A minimum-uncertainty 
solution of Eq. ([0j) would have Z = 0. 

Another interesting issue is the relationship between the ESS parameters a, ft, 
7o, 71, 5 and quantities determining the shape of the corresponding classical orbital 
ellipse, namely, the semimajor axis a, the angle r] it makes with the x axis, and the 
eccentricity e. The Ehrenfest theorem requires that the mean values of the position 
(r) and the momentum (p) evolve according to the equations d(f)/dt = (p) and 
d(p)/dt = -(W(f)). This means that the mean position moves along the classical 
trajectory if (VV(f)) = W((f)). However, this condition is not satisfied for the 
Coulomb problem. The connection between the ESS parameters and the shape of 
the corresponding classical orbit is therefore difficult to establish analytically and is 
partly a matter of definition. 
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If the quantum motion did follow the classical trajectory, knowledge of the initial 
position and momentum would suffice to determine it. One useful approximate re- 
lationship between the ESS parameters and those of the classical orbit in cartesian 
coordinates is therefore given by taking ESS expectations of suitable quantities as 
the initial data for a classical trajectory and thereby establishing the desired relation. 
This procedure has the advantage of being analytical and providing insight about the 
meaning of the ESS parameters. 

For present purposes, we choose to express the ESS parameters as functions of a, 
rj, e, and the widths Ar and 5 of the ESS, which have no direct classical counterpart. 
We find 

Iq(25) a ( e 2 sin 2 rj \ 

A (2d) 2(Ar) 2 \ 1 — ecosrjJ 

a = 2 7o 2 (Ar) 2 -l , (47) 



(2« + l) 2 , n . (I {28)' 



2a + 1 . (h{25)' 



3 



-» = -m^rr^Ws)) • (49) 

Note that there are two solutions to the equation for (3, reflecting the two possible 
directions the trajectory is followed. For the examples considered above, rj = by 
assumption, corresponding to an ESS located at the outer apsidal point. 

As an example of the use of these expressions, we can address the general issue 
of whether Eq. ( fH|) is minimized by an ESS for any shape of classical orbit. For 
definiteness, choose Ar and AL to match the example in Fig. 2, and take rj = 0. 
Then, the expressions ( f4"9"D can be used to obtain the quantity Z of Eq. fllPp as a 
function of the semimajor axis a and eccentricity e. Figure 4 shows a plot of the 
resulting lengthy expression for Z, which provides a normalized measure of closeness 
to minimum uncertainty in Eq. (|44]). The value Z = corresponds to minimum 
uncertainty. Figure 4 shows that Z > for the full range of associated classical orbits 
displayed in the graph. This demonstrates more generally that the ESS are different 
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from the superposition of 0(3) coherent states presented in ref. ||. 



V. INCORPORATION OF QUANTUM DEFECTS 

Most experiments studying the properties of Rydberg wave packets are performed 
using alkali-metal atoms. These have quantum defects causing energy-level shifts 
away from the hydrogenic values. It is therefore useful to have a relatively simple 
and analytical theory for alkali-metal atoms that has eigenstates with the attractive 
properties of hydrogenic eigenf unctions such as completeness and orthogonality, but 
that generates exact quantum-defect energies. A model satisfying these criteria, called 
supersymmetry-based quantum-defect theory (SQDT), is known to exist [Q. In this 
section, we outline the planar analog of SQDT and sketch the construction of the 
associated ESS. These may be of use in modeling packets produced in alkali-metal 
atoms. 

The physical quantum defects 5(n, I) are empirical parameters depending in gen- 
eral on n and /. We define the planar analogues S(n, |/|) as functions of the absolute 
value so that a planar Rydberg series can be established having features in com- 
mon with the standard case. As usual, the quantum defects for large n approach 
asymptotic values 8(\l\) that are independent of n. 

The planar Rydberg series has the form 

E "- = W^W • (50) 

where n* = n — 5(\l\). The quantum defects partially lift the degeneracy among 
states having different values of |/|. The primary objective of planar SQDT is to 
obtain analytical eigenfunctions with ([5(]) as eigenenergies. This is accomplished by 
defining an effective radial potential 

1 \l*\ 2 -\l\ 2 , , 

where \l*\ = \l\ —S(\l\) + and is an integer that depends on The radial 
eigenstates R n *i*(r) of this potential have the form of R n i in Eq. fl39|) but with n and 
\l\ replaced by n* and \l*\, respectively. For certain values of the integers /(|/|) and 
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in the limit of vanishing quantum defects, these solutions admit a supersymmetric 
extension. 

The set of planar SQDT eigenstates R n *i*(r)Yi(<f)), where Y\ is an eigenstate of 
angular momentum, is complete and orthonormal. We may choose S(\l\) to match 
the asymptotic quantum defects 5(1), with / > 0, for the physical alkali-metal atom. 
For example, for lithium s and p states, 5(0) ~ 0.40 and 5(1) ~ 0.05. With these 
values, the planar SQDT eigenstates reproduce the planar Rydberg series for lithium. 

An interesting issue is the ESS construction in the context of this model. At the 
classical level, the presence of the modified 1/r 2 term in the effective potential ( |51~D 
means the Runge-Lenz vector A is no longer conserved. The orbit precesses at a 
rate determined by the quantum defect and is no longer a keplerian ellipse. At the 
quantum level, the commutator of the Runge-Lenz vector with the hamiltonian is 
nonzero, and the SO (4) symmetry is broken. 

Despite these differences, the RSS construction can be performed in the presence 



of quantum defects For the present case, the oscillator variables R and P must 
be modified to 

R = -~tL . P = -Udr + Pl , (52) 



r |/*| 2 ' /V 2r 

where / = |/*|/|/|. These operators obey the uncertainty relation 

AJ2AP>i(i> , (53) 
2/ r z 

which also depends on the quantum defect. The wave function ip(r) minimizing this 
relation has the same functional form as that in Eq. (|5|). 

An ESS \l/(r, (ft) can be formed as a product of an RSS ip(r) and a CSS The 
form of the solution is like that in Eq. fl2"E|). However, some of the parameters must 
be chosen differently. We choose (p r ) = and (L) = I, and we specify AL as before. 
However, (H) = E n * and (r) = r* ut differ from the previous case. Here, E n * is the 
energy expectation and r* ut is the outer apsidal point for a superposition of states 
with quantum-defect eigenenergies. This means the ESS is initially located at the 
outer apsidal point of a precessing ellipse. The time evolution and revival structure 
depend on the quantum defects and can be studied by expanding \l/(r, 0, t) in the 
complete set of SQDT eigenstates. 
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In the above, a prescription for calculating (H) is needed since the SQDT potential 
depends on |Z| and |/*|. We can take advantage of completeness to expand the initial 
ESS as a superposition of SQDT eigenstates: 

^(r,0) = ]Tg n ^,(r)YK0) , (54) 

n,l 

where the expansion coefficients c n u which depend on the ESS parameters, can be 
determined by inversion. The expectation value for the hamiltonian can then be 
specified as 

(H) = J2\cm\ 2 E n ,=E n , . (55) 

n,l 

VI. SUMMARY 

In this paper, we have found analytical solutions to the planar Coulomb problem, 
called elliptical squeezed states, that minimize coordinate-momentum uncertainty 
relations and that move along classical keplerian orbits. The paper also provides an 
extension of the analysis to the case where quantum defects are present, which is 
of experimental importance but theoretically difficult to treat in other approaches 
because quantum defects break the 0(4) symmetry. 

The ESS provide an analytical tool for studying the quantum-classical corre- 
spondence in the Coulomb problem. They may also be used to describe minimum- 
uncertainty Rydberg wave packets excited by short laser pulses with external fields 
present. To match an ESS to a Rydberg wave packet, we initialize the ESS at the 
outer apsidal point of the orbit. The five ESS parameters are given in terms of the ex- 
pectation values of the radial position, the radial momentum, the energy, the angular 
momentum, and the spread in the angular momentum. 

We have studied the time evolution of the ESS for examples with different semi- 
major axes and eccentricities. The wave packets move along an elliptical trajectory 
with the classical keplerian orbital period. The squeezing of the ESS is evident as the 
width of the wave packet oscillates during the motion. The ESS maintain their form 
for several orbital cycles before decoherence. 
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FIGURE CAPTIONS 



Fig. 1: Sample CSS |x(0)| 2 (arbitrary units) plotted as a function of 4> in radians. 
All three cases shown are centered on <fi = with f3 = (L) = 30. The angular- 
coordinate widths of the CSS depend on S and have the values AL = 0.5 for 
5 ~ 0.804 (solid line), AL = 1.5 for 5 ~ 4.757 (dotted line), and AL = 2.5 for 
5 ~ 12.753 (thick solid line). 

Fig. 2: The modulus of an ESS r|\&(r, 0, t)| 2 (arbitrary units) plotted as a function 
of r and at the times: (a) t = 0, (b) i = |T c i, (c) t = |T d , (d) t = |T cl , (e) 
t = T cl . 

Fig. 3: Comparison of the modulus squared (arbitrary units) of two ESS with different 
eccentricities at the time t = \T C \. Both wave packets have n = 45 and AL = 
2.5. and are viewed from a point on the positive looking in toward the 

origin. For Figure (a), I — 30, while for Figure (b), / = 40. 

Fig. 4: A plot of the quantity Z = (AA X AA y - \(H L)\) / \(H L)\ as a function of 
the semimajor axis a in atomic units and the eccentricity e of the associated 
classical orbits. 
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